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bstract — In this paper, we investigate optimal coding strate- 
gies for a class of linear deterministic relay networks. The 
network under study is a relay network, with one source, 
one destination, and two relay nodes. Additionally, there is a 
disturbing source of signals that causes interference with the 
information signals received by the relay nodes. Our model 
captures the effect of the interference of message signals and 
disturbing signals on a single relay network, or the interference of 
signals from multiple relay networks with each other in the linear 
deterministic framework. For several ranges of the network 
parameters we find upper bounds on the maximum achievable 
source-destination rate in the presense of the disturbing node 
and in each case we find an optimal coding scheme that achieves 
the upper bound. 

I. Introduction 

Finding the capacity of relay networks is a basic open prob- 
lem in network information theory (TJ. One way to approxi- 
mate the function of relay networks is through deterministic 
models that was first introduced by Aref (2). Avestimehr, 
Diggavi and Tse have recently Q, PJ introduced a linear 
deterministic model for wireless relay networks that treats 
noise as a deterministic thresholding function and the inter- 
ference of signals as a linear transformation over a finite field. 
They have successfully applied this model to several relay 
networks and have shown that the capacity of the deterministic 
model is within a constant gap from the corresponding wireless 
network. Furthermore, they have shown that a max-flow min- 
cut result holds for the capacity of the relay deterministic 
network in the case of a single multicast session. While the 
capacity achieving scheme in ll3l. PI is a random coding 
over long blocks of signals, recent works 0, J6), Q, I0 
have devised low complexity and deterministic schemes that 
achieves the maximum capacity in the case of a unicast 
session. 

In the case of multiple messages, Mohajer et al., J51 have 
considered two unicast sessions on a relay network with two 
relays, two sources, and two destinations. They study two 
special cases of this setting, namely ZS and ZZ channels and 
give a full characterization of the capacity region along with 
capacity acheiving schemes in each case. 

In this paper, we are also interested in the deterministic 
relay network with interference at each node. The network 
here consists of a source node, a destination node, and two 
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relay nodes. Additionally there is a disturbing node that sends 
signals to the two relays and causes interference with message 
signals. Our goal is to characterize the capacity region from 
source to destination in the presense of the disturbing node. 
Our model differs from the model in (9) in two ways. First 
they consider a two dimensional capacity region in which 
each source node tryes to send messages to its designated 
destination. But in our model, we only try to find the one 
dimensional capacity from one source to the corresponding 
destination. On the other hand, while in both the ZS and ZZ 
channels the interference of the two messages only occur in 
one relay node, in our model we assume that the disturbing 
signals have interference with both relays which makes it more 
difficult to handle. We find the capacity and offer capacity 
achieving schemes that have low complexities for this relaying 
problem. 

The organization of the material in this paper is as follows. 
In Section II we describe the deterministic model of relay 
networks and the model of the network that we study in 
this paper. Furthermore, we discuss the general linear coding 
and decoding strategies and the achievable rate of the linear 
schemes. In Section III we find optimal linear coding schemes 
and their corresponding achievable rates. 

II. Problem Setting 

First we briefly state the linear deterministic model of relay 
networks from Q, RJ, then we introduce our considered 
network model. 

A. Deterministic model 

Consider a directed graph N(V, E) where V denotes the 
set of nodes in the network including source, relays and 
destination and E is the set of edges. Communication from 
node i to j has a nonnegative gain n,j associated with it. 
This number models the channel gain in the corresponding 
Gaussian setting. Each node i transmits a vector xi £ F% 
and receives a vector 6 F2 where q = rnaxj^ -.(naj-)). 
The received signal at each node is a deterministic function 
of transmitted signals at the other nodes with the following 
input-output relation: 
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Figure 1. Diamond network with a disturbing node 

where Q is the q x q shift matrix given by 
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In this paper, we are interested in linear coding schemes where 
at every node j the transmitted signal is a linear function of 
the received signal 

x j = G jVj- (2) 

The goal is to design coding functions Gi such that the 
destination node receives enough information for decoding the 
message sent by the source. 

B. Diamond network with a disturbing node 

We consider a diamond relay network with a disturbing node 
which is depicted in Figure [T] Here S and D are respectively 
the source and the destination nodes and M is the source of 
disturbing signals. Also nodes A and B are the relay nodes. 
For simplicity of our analysis we assume that the links from 
M to the two relays have the same gain m that is also realistic 
in the situation where the distance between the nodes in the 
diamond wireless network is relatively small compared to the 
distance from node M which might be operating in a different 
network. As before, we denote the transmitted signal from 
node i by x$ and the received signal by x/i. 

Let Ga and Gb be the coding matrices at nodes A and B 
respectively. The signal y B received by destination is a linear 
combination of the transmitted signals xs and xm- Let 

Vd = G s x s + G m xm- 
where by application of ([TJ we have 

G s = Q q - ni G A Q q - ni + Q q - ni G B Q q - n \ (3) 
G M = Q q - n3 G A Q q - m + Q q - n *G B Q q - m . (4) 

For every choice of Ga and Gb we let = R(Ga,Gb) 
denote the rate of transmission of information from S to D. 
Also we let C = maxc A .G B R(Ga,Gb) denote the capacity 
of transmission from S to D. In the following, we explicitly 
find R(Ga, Gb) and in the next section we find the capacity 
C in terms of the network parameters. 



In our analysis, we usually work with the range and di- 
mension of matrices. For matrix H we let range (H) denote 
the linear span of the columns of H. Also for a subspace 
S C ¥ 2 we let dim(iS) denote its dimension. Obviousely 
dim(range(i/)) = rank(i7). For two matrices Hi, H2 we also 
use the shorthand rank(ii! n H 2 ) to denote dim (range (Hi) n 
range (H 2 )). 

Theorem 1. For any choice of Ga and Gb we have 

R(G A ,G B ) = rank(G s ) - rank(G s n G M ) 

where Gs and Gm are defined in and Q. 

For two subspaces S\ and S 2 of F|, let Si + S 2 = 
{si + S2 : Si G Si, S2 G 52}. To prove the Theorem[T]we first 
prove the following lemma. 

Lemma 1. If Si and S 2 are two subspaces of ¥ 2 then for 
any s G S = S\ + £2 there exists a unique pair Si G Si and 
s 2 G S 2 with si + s 2 = s if and only if Si D S 2 = {0} . 

Proof: First suppose that Si <~) S 2 ^ {0} hence, there 
exists at last t 7^ such that t G Si , t £ S 2 . Then if s = si+s 2 
with si G Si and s 2 G S 2 , we also have si + t G Si and 
s 2 — t G S 2 and (si + t) + (s 2 — t) = s. Therefore the 
condition is necessary. To prove the sufficiency, we notice that 
S is a subspace of ¥ 2 and if Si S 2 = {0} , we can form 
a basis for S by union of a basis of Si and a basis of 6>2. 
Now if s = si + s 2 = t\ + t 2 and si ^ t\, s% ^ t 2 we can 
form two different expansions for s in the basis of S by using 
either the expansion of s^nd s 2 in the bases of Si and 52 or 
the expansions of ti and t 2 in the bases of Si and 6>2- This 
contradicts the fact that each s G S has a unique expansion in 
the basis of S. ■ 
Next we prove Theorem [T] 

Proof: For a linear scheme we choose a subspace X C F$j 
as the set of our codewords x$- Let S = {G$xs '■ x s G ^} ■ 
For a successful decoding X has to satisfy two properties; 
first, each xs G X should be mapped into a unique vector in 
S. This implies that dim(5) = dim(A'). Second, for every 
x s G X, and every x M G F^j, y D = G s x s + G M x M 
corresponds to the unique pair of Gsxs and Gm x m- These 
two conditions guarantee that each yD corresponds to a unique 
codeword xs- Then, the maximum dimension of X that 
satisfies the two conditions is the rate R of the code. Since 
Gm%m can be any vector in range(GA/ ), to satisfy the second 
condition, by Lemma [T] range(Gjif) H S = {0} . S has the 
maximum dimension when it is the largest subspace of the 
set (range(Gs) — range(GA/)) U {0} which has a dimension 
of rank(Gs) — rank(Gs n Gm). Next we choose X to be the 
subspace of ¥ 2 with dimension of rank(Gg) — rank(GsnGM) 
such that Gs maps it to the set S. Notice that since S is a 
subspace of range^g) we can always find such X. Therefore 
R = dim{X) = dim(S) = rank(G s ) - rank(G s n G M ). ■ 

III. Capacity of the Network 

In this section, we find the linear capacity of the network, 
that is the maximum achievable rate by a linear coding scheme. 



In several steps of our capacity calculation we will use the 
following useful lemma: 

Lemma 2. Let F mxn and G rnxn be two matrices that are the 
same in at least their first to — a rows. Then 

rank(.F) — rank(F n G) < min(n, a). 
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rank(F) - rank(F D G) = rank( [ F G ] ) - rank(G). 
We have 

rank([ F G } ) = rank( ^ ^ 
rank([A]) + rank( [ B D ] ) < rank(A) + a. 

A 
D 



Also, rank(G) 



rank( 



) > rank(A). Therefore 



rank( [ F G ] ) - rank(G) < a. Also, rank(F) < n and 
thus rank(F) - rank(F n G) < n. ■ 

Next we derive the linear capacity of the network in Figure 
[T] By symmetry, we only derive the capacity for ri\ >ri2- 

1) ni > ri2, U3 > ri4 and m < n\ : Let A be a q x q 
matrix. For two integers ni and n r consider a partition of the 
elements of A as follows 
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Now let Ga be partitioned as follows 
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Also 

Q q - n3 G A Q q - m 

Q q ~ ni G B Q q - m 
Therefore, 



q - "3 

G M = n 3 -n 4 

"4 

Next we find an upper bound on any acheivable rate 
R = rank(Gs) — rank(Gs P| Gm) m terms of the network 
parameters and then find Ga and Gb that achieve the bound. 





R = rank(Gs) - rank(G s n G M ) 
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< ni — to + min(m, n 4 ), 

where (a) follows by Lemma [2] On the other hand, R can not 
exceed the value of any cut set in the absence of the noisy 
source, hence 

R < max(ri3, n 4 ) = n 3 

if we let k — min(m, n 4 ) then we have 

R < min(n! — m + k, n 3 ). (5) 

We let matrices Ga and Gb be of the following forms 
and show that they achieve the above upper bound. Let 
j = min(ni — m,n 3 — k) and for any t let it be the t x t 
identity matrix. 
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Then it is easy to verify that G5 is of the following form 
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where L is a lower triangular matrix with ones on its diagonal. 
This implies that rank(Gs) = j + k. Also Gm = qX q and 
hence a rate of R = j + k is achievable by this coding scheme. 
This matches our bound ([5). 

2) m > ri2, n 3 > ri4 and to > m: In this case Gs and 
Gm are as the following formats: 



Let us define 
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We have the following upper bound on any achievable rate 
R: 



R = rank(Gc 



< rank( 
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where (a) follows by Lemma [2] Let k = min(7ii, 714). We 
design Ga and Gs as follows: 
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= Ogxq and Gs is of the following form 
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where L is a lower triangular matrix with ones on its diagonal. 
Therefore rank(Gs) = k and R — k is achievable. This 
matches bound {6). 

3) rii > 712, 7^4 > ^3 ant/ m < n 2 : In this case for 
any choices of Ga and Gb,Gs and Gm are of the following 
forms: 
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The following bound holds for any rate R : 
R = rank(G s ) - rank(G s n G M ) 
< rank( 
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where (a) follows by Lemma [2] This bound, together with the 
cutset bound R < min(rti, 124, n 2 + n 3 ) in the absence of the 
node A/, results in the following bound 

R < min(ni,n4,ri2 + n 3 — m). (7) 

Let k = min(rt3,TO.). We design Ga and Gb as follows: 
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F : 

We know that: 

max (rank (i 71 )) = min(n4 — k, ni — to, 712 + n 3 — to — k). 
By this choice, Gm = gX9 and G5 is of the following form 



G s = 



where L is a lower triangular matrix with ones on its diagonal. 
Therefore 

rank(G s ) =k + mnk(F) 

= min(ri4 — k, rii — to, n 2 + n 3 — to — k) + k 
= min(?i4, ni — to + k, n 2 + n 3 — to) 
= min(n4, ni — to + n 3l ni — m + to, ^2 + n 3 — m) 
= min(ni , 714 , n 2 + n 3 — to) 

and i? = min(ni, 77,4, n2+n3 — to) is achievable. This matches 
our bound (|7j. 



4) ni > n 2 , «4 > n 3 and m > n 2 : In this case Gs and 
Gm are of the following forms 
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Let us define the following matrix 



Gm — 



Then we have the following bound on R: 

R = rank(Gs) - rank(G s n G M ) 
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where (a) follows by Lemma [2] Let k — min(ni,n3). We 
design Ga and Gb as follows: 
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It is easy to verify that Gm = Okxk and Gs is of the following 
form 

Tlx — k k q — ri\ 
q-k ( 
k { * L 



G s = 



where L is a lower triangular matrix with ones on its diagonal. 
Therefore R = rank(Gs) = k. This matches the bound (j8]l. 

5) ni — n 2 , m > rii : In this case it is easy to verify that 
for any choice of Ga and Gb, columns of Gs are a subset of 
columns of Gm- Therefore rank(Gg fl Gm) — rank(Gs) and 
hence R = rank(G s ) - rank(G s n G M ) = 0. 

6) ni — n 2 , m < n\: In this case for any choice of 
Ga and Gb, columns of Gm are a subset of columns of G5. 
Therefore rank(GsflGAf ) = rank(GM) 7 and R — rank(Gs) — 
rank(GM)- Since there are at most n-y — m columns in Gs 
that does not appear in Gm, rank(Gs) — rank(GM ) < n\ —m 
and therefore R < n\ — m. Also from the cutset bound R < 
max(?"i3, 71,4) we have 



Let k = min(ni — m, max(n 3 , TI4)). Consider two cases: 
• If n>3 > ^4 then we set Gb — Oqxq and let 
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Therefore R = rank(Gs) — rank(GM) = k, that achieves 
bound 

If TI4 > n 3 then we set Ga = qxq and let 
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Similar to the previous case, in this case Gm = Oqxq 
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R < min(ni — m, max(ri3, 714)). 



(9) 



Therefore R = rank(G5) — rank(GM) = k, that achieves 
bound |9]). 

Remark 1. It is easy to verify that for cases 1, 3 and 6 where 
m can be set to zero, the capacity of diamond network Q is 
achievable by our coding scheme. 
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